%%Valutazione frequenza cardiaca
%%RR = Tempo in millisecondi della distanza tra due picchi R
HR = 60./(RR./1000);
figure;
plot (HR);

%Media su primi 20 RR
numMed = 20;
HRIniziale = 0;
for (i=1:1:numMed)
    HRIniziale = HRIniziale + RR(i,1);
end
HRIniziale = HRIniziale/numMed

c = max(size(RR));
HRFinale = 0;
for (i=1:1:numMed)
    HRFinale = HRFinale + RR(c-i,1);
end
HRFinale = HRFinale/numMed

HRIniStr = ['HR Iniziale = ' num2str(HRIniziale)];
HRFinStr = ['HR Finale = ' num2str(HRFinale)];
title (['\fontsize{15} \color{red} \it      Hearth rate' '\newline \fontsize{10} \color{black}' HRIniStr '; ' HRFinStr]);
grid on;